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Abstract 

Participation  factors  are  nondimensional  scalars  that  measure  the  interaction  be¬ 
tween  the  modes  and  the  state  variables  of  a  linear  system.  Since  their  introduction 
by  Verghese,  Perez-Arriaga  and  Schweppe,  participation  factors  have  been  used  for 
analysis,  order  reduction  and  controller  design  in  a  variety  of  fields.  In  this  paper, 
participation  factors  are  revisited,  resulting  in  new  definitions.  The  aim  of  these  defi¬ 
nitions  is  to  achieve  a  conceptual  framework  that  doesn’t  hinge  on  any  particular  choice 
of  initial  condition.  The  initial  condition  is  modeled  as  an  uncertain  quantity,  which 
can  be  viewed  either  in  a  set-valued  or  a  probabilistic  setting.  If  the  initial  condition 
uncertainty  obeys  a  symmetry  condition,  the  new  definitions  are  found  to  reduce  to 
the  original  definition  of  participation  factors. 

Keywords:  participation  factors,  linear  systems,  modal  analysis,  stability,  dynamics,  prob¬ 
ability 

1  Introduction 

Since  its  introduction  by  Verghese,  Perez-Arriaga  and  Schweppe  [9], [5], [10],  Selective  Modal 
Analysis  (SMA)  has  become  a  popular  tool  for  system  analysis,  order  reduction  and  actuator 
placement.  In  particular,  this  tool  is  extensively  used  in  the  electric  power  systems  area  [4], 
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Participation  factors,  a  key  element  of  SMA,  provide  a  mechanism  for  assessing  the  level  of 
interaction  between  system  modes  and  system  state  variables. 

In  this  paper,  participation  factors  are  revisited,  resulting  in  new  definitions.  The  pur¬ 
pose  of  these  definitions  is  to  achieve  a  conceptual  framework  that  doesn’t  hinge  on  any 
particular  choice  of  initial  condition.  The  initial  condition  is  modeled  as  an  uncertain  quan¬ 
tity,  which  can  be  viewed  either  in  a  set-valued  or  a  probabilistic  setting.  If  the  initial 
condition  uncertainty  obeys  a  symmetry  condition,  the  new  definitions  are  found  to  reduce 
to  the  original  definition  of  [9], [5], [10].  Since  the  initial  condition  is  viewed  as  uncertain, 
the  definitions  model  its  effect  on  the  level  of  mode-state  interaction  in  an  average  sense. 
(For  simplicity,  the  term  “state”  is  used  here  interchangeably  with  “state  variable.”)  This 
work  provides  a  framework  for  a  deeper  appreciation  of  participation  factors,  as  well  as  an 
opening  to  further  useful  generalizations  of  the  concept  in  various  directions.  For  instance, 
generalizations  to  other  system  types  can  be  pursued.  Also,  definitions  along  the  same  lines 
could  facilitate  analytical  treatment  of  the  relation  between  control  inputs  and  states/mode 
interaction,  which  in  turn  would  have  implications  for  actuator  and  sensor  placement. 

Fowllowing  Verghese,  Perez- Arriaga  and  Schweppe  [9], [5], [10],  participation  factors  are 
considered  in  two  basic  senses.  In  the  first  sense,  a  participation  factor  measures  the  relative 
contribution  of  a  mode  to  a  state.  In  the  second,  a  participation  factor  measures  the  relative 
contribution  of  a  state  to  a  mode.  It  isn’t  clear  at  the  outset  that  these  two  senses  should  lead 
to  identical  formulas  for  participation  factors.  However,  the  precise  definitions  in  [9], [5], [10] 
for  these  two  senses  of  participation  factors  did  indeed  result  in  identical  mathematical 
expressions.  The  same  conclusion  is  found  to  apply  in  the  present  paper,  under  assumptions 
valid  for  a  large  class  of  problems. 

It  should  be  noted  that  there  have  been  other  interpretations  of  participation  factors 
since  the  original  work  of  [9], [5], [10].  For  example,  participation  factors  are  often  viewed  as 
sensitivities  of  eigenvalues  to  changes  in  the  diagonal  entries  of  the  state  dynamics  matrix 
(see,  e.g.,  [8]).  Interpretations  in  terms  of  modal  energies  are  also  common  (see,  e.g.,  [1]). 
Still  another  interpretation  relates  to  eigenvalue  mobility  under  state  feedback  [6]. 

The  remainder  of  the  paper  is  organized  as  follows.  In  Section  2,  the  original  definition  of 
participation  factors  [9], [5], [10]  is  recalled,  and  motivation  for  the  work  of  this  paper  is  given. 
In  Section  3,  the  new  definitions  of  this  paper  that  address  participation  of  modes  in  states 
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are  given.  The  relationship  with  the  original  definition  is  shown.  Section  4  contains  the  new 
definitions  that  address  participation  of  states  in  modes  along  with  their  relationship  to  the 
original  definition  in  [9], [5], [10].  Concluding  remarks  are  given  in  Section  5. 

2  Background  and  Motivation 

Participation  factors  were  introduced  by  Verghese,  Perez- Arriaga  and  Schweppe  [9], [5], [10] 
as  a  means  for  ranking  the  relative  interactions  between  system  modes  and  system  states. 
The  concept  is  one  element  of  the  Selective  Modal  Analysis  (SMA)  approach  introduced  by 
these  authors,  and  its  first  applications  were  in  the  field  of  electric  power  systems.  This 
section  recalls  the  original  definition  of  participation  factors. 

Consider  a  general  continuous-time  linear  time-invariant  system 

x(t)  =  Ax(t ),  (1) 

where  x  G  Rn,  and  A  is  a  real  n  x  n  matrix.  Suppose  for  simplicity  that  A  has  a  set 
of  n  distinct  eigenvalues  (Ai,  A2, . . . ,  An).  Typically  the  evolution  of  each  state  variable  is 
influenced  by  all  the  eigenvalues  (Ai,  A2, . . . ,  An).  However,  it  is  often  desirable  to  quantify 
the  participation  of  a  particular  mode  (i.e.,  eigenmode)  in  a  state  variable.  If  the  state  vector 
is  composed  of  physically  meaningful  variables,  such  a  quantification  leads  to  conclusions 
regarding  the  influence  of  system  modes  on  physical  components. 

It  is  at  first  tempting  to  base  the  association  of  modes  with  state  variables  on  the  magni¬ 
tudes  of  the  entries  in  the  right  eigenvector  associated  with  a  mode.  Let  (r1,r2, . . . ,  rn )  be 
right  eigenvectors  of  the  matrix  A  associated  with  the  eigenvalues  (A1;  A2,  •  •  • ,  An),  respec¬ 
tively.  Using  this  criterion,  one  would  say  that  the  mode  associated  with  A*  is  significantly 
involved  in  the  state  Xk  if  rlk  is  large.  As  pointed  out  in  [9], [5], [10],  this  approach  has  two 
main  disadvantages.  The  first  is  that  it  requires  a  complete  spectral  analysis  of  the  system, 
and  is  thus  computationally  expensive.  The  second,  which  is  the  more  serious  flaw,  is  that 
the  numerical  values  of  the  entries  of  the  eigenvectors  depend  on  the  choice  of  units  for  the 
corresponding  state  variables.  This  renders  the  criterion  unreliable  in  providing  a  measure 
of  the  contribution  of  modes  to  state  variables.  This  is  true  even  if  the  variables  are  similar 
physically  and  are  measured  in  the  same  units. 
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The  foregoing  was  the  motivation  for  the  approach  taken  in  SMA  to  quantify  the  rel¬ 
ative  contributions  of  modes  to  state  variables.  In  SMA  the  entries  of  both  the  right  and 
left  eigenvectors  are  utilized  to  calculate  “participation  factors”  that  measure  the  level  of 
participation  of  modes  in  states  and  the  level  of  participation  of  states  in  modes.  The  par¬ 
ticipation  factors  defined  in  SMA  are  dimensionless  quantities  that  are  independent  of  the 
units  in  which  state  variables  are  measured. 

Next,  a  brief  summary  of  the  original  definition  of  participation  factors  in  the  sense 
of  participation  of  modes  in  states  is  given.  As  noted  earlier,  the  original  definition  also 
encompasses  participation  of  states  in  modes. 

Let  (l1,  l2, . . . ,  ln )  denote  the  left  (row)  eigenvectors  of  the  matrix  A  associated  with  the 
eigenvalues  (Ai,  A2, . . . ,  An),  respectively.  The  right  and  left  eigenvectors  are  taken  to  satisfy 
the  normalization  [3,  p.  154] 


llr3  =  8ij:  1  <i,j<n, 


where  Sij  is  the  Kronecker  delta  symbol: 


Sij  = 


1  if  i  =  j, 
0  if  i  ^  j. 


(2) 


In  the  remainder  of  the  paper,  the  i-th  mode  refers  to  the  mode  associated  with  Xt.  indeed 
specifically  to  the  eXit  terms  in  the  state  trajectory  (and  not,  for  simplicity,  to  the  conjugates 
of  these  terms  in  case  of  a  complex  eigenvalue;  see  Remark  2  below). 

In  [9], [5], [10],  the  participation  factor  of  the  i-th  mode  in  the  k- th  state  is  defined  to 

be  the  complex  number 


Pki  ■=  4rfe-  (3) 

The  motivation  for  this  definition  given  in  [9], [5], [10]  is  as  follows. 

The  solution  of  the  dynamic  system  equation  (1)  satisfying  the  initial  condition  x(0)  =  x° 

is 


x(t)  =  eMx°.  (4) 

Since  the  eigenvalues  of  A  are  assumed  distinct,  A  is  similar  to  a  diagonal  matrix.  Using 
this,  Eq.  (4)  gives 

n 

x(t)  =  ^(fx°)eAiV .  (5) 

i= 1 
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This  formula  is  alternatively  (and  conveniently)  viewed  as  the  superposition  of  solutions 
corresponding  to  initial  conditions  along  the  eigenvectors,  where  each  such  initial  condition 
is  the  projection  of  x°  along  the  eigenvector. 

Now  suppose  the  initial  condition  x°  is  ek .  the  unit  vector  along  the  k- th  coordinate  axis. 
Then  the  evolution  of  the  /c-th  state  becomes 

n 

xk(t)  =  J2(lkrl)eXit 

i= 1 
n 

=  Y,PkieXit.  (6) 

i= 1 

Eq.  (6)  indicates  that  pki  can  be  viewed  as  the  relative  participation  of  the  i-th  mode  in  the 
k- th  state  at  t  —  0. 

Remark  1.  The  quantities  defined  in  Eq.  (3)  are  also  used  in  [9], [5], [10]  to  measure  rela¬ 
tive  participation  of  states  in  modes.  The  motivation  for  this  interpretation  given  in  [5]  is 
somewhat  similar  to  the  above  and  need  not  be  reproduced  here.  ■ 

As  noted  above,  any  notion  of  participation  factors  is  useful  only  if  it  results  in  quantities 
that  do  not  depend  on  the  units  in  which  state  variables  are  measured.  In  this  remark,  we 
record  for  later  use  the  effect  of  changes  in  units  on  right  and  left  eigenvectors.  A  general 
change  in  units  is  represented  by  a  transformation  Xj  =  ajXj,  j  —  1, ...  ,n,  where  the  otj  are 
positive  constants.  Let  rl  and  ll  be  right  and  left  eigenvectors,  respectively,  corresponding 
to  eigenvalue  A*,  in  the  original  units.  It  is  straightforward  to  show  that  these  eigenvectors 
become,  in  the  new  units,  fl  and  ll  with  components  r*  =  ct,r*  and  /*  =  aj 1  /*• .  In  particular, 
it  is  clear  that  the  quantities  defined  in  (3)  above  are  independent  of  units. 

The  approach  of  the  present  paper  to  defining  participation  factors  builds  on  the  original 
work  in  [9], [5], [10].  The  focus  is  on  extending  the  participation  factors  concept  of  [9], [5], [10] 
to  explicitly  incorporate  the  effect  of  uncertainty  in  the  initial  condition  x°. 

The  reconsideration  of  modal  participation  in  this  light  is  motivated  by  the  following 
simple  observations  for  the  linear  system  (1)  and  its  solution  (5).  If  the  initial  condition 
x°  lies  along  the  z-tli  eigenvector,  then  the  only  mode  that  participates  in  the  evolution  of 
any  state  is  the  i-th  mode  (i.e.,  the  eXit  mode).  On  the  other  hand,  if  the  initial  condition 
lies  along  the  /c-th  coordinate  axis,  then  the  evolution  of  the  /c-th  state  involves  all  system 
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modes  according  to  Eq.  (6)  above.  Clearly,  then,  the  degree  to  which  a  mode  participates 
in  a  state  depends  on  the  initial  condition.  Similar  considerations  can  be  given  regarding 
participation  of  states  in  modes. 

In  the  next  section,  a  new  approach  to  defining  participation  factors  is  given  for  the  case 
of  participation  of  modes  in  states.  Section  4  gives  the  corresponding  results  for  participation 
of  states  in  modes.  Both  sections  rely  on  explicitly  incorporating  the  effect  of  uncertainty  in 
the  system  initial  condition  x°  through  an  averaging  operation. 

3  Participation  Factors:  Modes  in  States 

The  linear  system  (1)  usually  represents  the  small  perturbation  dynamics  of  a  nonlinear 
system  in  the  neighborhood  of  an  equilibrium.  The  initial  condition  for  such  a  perturbation 
is  naturally  viewed  as  being  an  uncertain  vector  of  small  norm. 

There  are  several  ways  in  which  the  foregoing  comments  can  be  implemented  to  result  in 
definitions  of  participation  factors.  For  example,  one  can  suppose  that  the  initial  condition 
is  known  to  lie  in  a  specified  symmetric  set  in  Htn  and  calculate  “average  participations”  of 
modes  in  state  variables  by  averaging  over  the  set.  This  is  done  in  Subsection  3.1.  If  the 
initial  uncertainty  set  is  not  symmetric  according  to  Definition  1  below,  then  the  set-theoretic 
approach  of  Subsection  3.1  does  not  lead  to  a  useful  notion  of  participation  factors.  In  some 
applications,  such  as  ecological  modeling,  it  is  natural  to  restrict  the  initial  condition  to  lie 
in  a  sector  and  not  in  a  full  neighborhood  of  the  equilibrium.  A  formalism  that  allows  such 
settings  is  given  in  Subsection  3.2.  The  approach  of  Subsection  3.2  entails  modeling  the 
initial  condition  uncertainty  probabilistically.  This  technique  is  used  again  in  Section  4  to 
streamline  the  analysis  of  participation  of  states  in  modes. 

3.1  Set- Valued  Initial  Condition  Uncertainty 

In  this  subsection,  the  initial  condition  x°  is  taken  to  he  in  a  connected  set  S  containing  the 
origin: 

x°  e  5.  (7) 

In  fact,  the  case  of  greatest  interest  is  when  S  =  ]Rn. 

Sets  S  that  are  symmetric  in  the  sense  of  the  next  definition  are  of  particular  significance. 
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Definition  1.  The  set  S  is  symmetric  with  respect  to  each  of  the  hyperplanes  Xk  =  0, 
k  =  1, . . . ,  n.  That  is,  for  any  k  G  {1, . . . ,  n}  and  z  =  (zi, . . . ,  Zk, . . . ,  zn)  G  Rn,  z  G  S 
implies  that  (zi, ,  —zk,  . . . ,  zn)  G  S.  ■ 

Although  probabilistic  arguments  aren’t  used  until  the  next  subsection,  intuitively,  Def¬ 
inition  2  below  is  consistent  with  assuming  a  uniform  density  for  the  initial  condition  and 
computing  the  expectation  of  the  relative  contribution  of  the  i-th  mode  in  the  k- th  state  at 
time  t  =  0.  Note  that  the  average  contribution  at  time  t  =  0  of  the  i-th  mode  to  state  Xk 
vanishes  and  so  is  not  useful  as  a  notion  of  participation  factor: 


avg  ( llx°)rlk  =  0.  (8) 

i°G<S 


Definition  2.  The  participation  factor  for  the  mode  associated  with  A*  in  state  Xk  with 
respect  to  an  uncertainty  set  S  symmetric  according  to  Definition  1  is 


Phi 


avg 
x°  G  S 


i„0) 


(. llx 


•A;  7. 


(9) 


whenever  this  quantity  exists.  Here,  “avgxog<s”  is  an  operator  that  computes  the  average 
of  a  function  over  the  set  S.  (In  computing  the  implied  multidimensional  volume  integral, 
the  argument  function  is  undefined  for  x°k  =  0,  and  the  Cauchy  principal  value  [7]  of  the 
integral  is  to  be  used.  For  example,  if  the  uncertainty  set  is  Rn,  then  a  limit  is  taken  over  a 
symmetric  set  S  that  tends  to  Rn .  The  limit  is  taken  in  a  symmetric  way  as  the  states  tend 
to  oo  at  the  upper  integration  limits  and  to  0  at  the  lower  integration  limits.)  ■ 

This  quantity  measures  the  average  relative  contribution  at  time  t  =  0  of  the  i-th  mode 
to  state  Xk-  In  the  definition,  the  i-th  mode  is  interpreted  as  the  eXit  term  in  (5).  Also,  the 
denominator  on  the  right  side  of  (9)  is  simply  the  sum  of  the  contributions  from  all  modes 
to  Xk(t)  at  t  =  0: 

n 

x°k  =  J2(l3x°)r{.  (10) 

3  = 1 


Remark  2.  When  A,  is  complex  its  conjugate  is  also  an  eigenvalue,  and  one  naturally  views 
the  mode  as  consisting  of  both  the  eXit  term  and  the  eXit  term  in  (5).  Because  of  this,  it 
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would  also  be  reasonable  to  define  the  participation  factor  for  complex  modes  as  the  twice 
the  real  part  of  the  quantity  defined  in  (9),  i.e.,  2R e(llkrlk).  This  is  not  done  here  for  two 
reasons.  The  first  is  to  simplify  the  presentation.  The  second,  which  is  the  more  important  of 
the  two,  is  to  ensure  that  the  new  definitions  not  differ  from  the  original  one  merely  because 
of  a  convention  such  as  whether  a  mode  includes  just  one  exponential  term  or  also  includes 
its  conjugate.  ■ 

Next,  the  expression  on  the  right  side  of  Eq.  (9)  is  simplified  for  the  case  of  a  symmetric 
set  S  according  to  Definition  1.  Let 

Vol(S)  :=  f  dx°  (11) 

Jx°es 

denote  the  volume  of  the  set  S.  If  the  set  S  has  infinite  volume,  then  the  construction  below 
is  performed  for  a  finite  symmetric  subset,  and  then  a  limit  is  taken  as  discussed  in  Definition 
2.  From  Eq.  (9),  pki  is  given  by 


Pu  =  avg  1W1+  avg  £ 
x°  G  S  Xk  x°  G  5  7  =  1 

3  +  k 


xi 


=  Hrl  + 


k1  k  1  /  ™0 

*0g<S  .  xk 

3  =  1 


dxu/Vol(5) 


3  ~f~  k 


n  r  x3 

=  44+  E  44  Jx0es  ^dx°/voi(5) 

3  =  1 
3  +  k 


llrl 
Lk'  k- 


(12) 


The  last  step  follows  from  the  observation  that,  because  S  is  symmetric  according  to  Defi¬ 
nition  1, 


x 


lx°es  xk 


^dx°  =  0 


(13) 


for  any  j  ^  k,  where  the  integral  is  interpreted  in  the  sense  of  the  Cauchy  principal  value. 
Summarizing  this  result,  it  has  been  found  that  for  a  set-valued  uncertainty  for  the  initial 


condition  with  a  symmetric  uncertainty  set  S,  Definition  2  for  pt,  reduces  to 

Pki  ■=  lkrk-  (14) 

This  result  for  p^  agrees  with  the  original  definition  of  Verghese,  Perez-Arriaga  and 
Schweppe  [9], [5], [10].  However,  Definition  2  is  based  on  the  symmetry  condition  as  given 
in  Definition  1.  Without  this  symmetry  condition,  the  integral  implied  in  the  definition  (9) 
becomes  cumbersome  to  evaluate,  and  its  Cauchy  principal  value  may  fail  to  exist. 

The  symmetry  assumption  on  S  is  reasonable  for  typical  engineering  system  models. 
However,  there  are  important  examples  of  system  models  for  which  the  state  vector  is  re¬ 
stricted.  For  example,  in  population  dynamics  and  chemical  process  dynamics,  state  vari¬ 
ables  are  often  restricted  to  be  nonnegative.  If  the  origin  is  an  equilibrium  for  such  a  system, 
then  initial  conditions  near  the  origin  must  be  restricted  to  have  nonnegative  components. 
In  part  to  accommodate  such  settings,  a  probabilistic  framework  to  studying  participation 
of  modes  in  states  is  pursued  in  the  next  subsection.  This  probabilistic  framework  will  be 
found  in  Section  4  to  greatly  facilitate  the  study  of  participation  of  states  in  modes.  Indeed, 
the  set- valued  approach  was  found  to  be  cumbersome  for  the  study  of  participation  of  states 
in  modes,  and  will  not  be  discussed  in  that  context. 


3.2  Random  Initial  Condition  Uncertainty 


Next,  a  probabilistic  definition  of  participation  factors  for  participation  of  modes  in  states 
is  given.  This  is  achieved  by  letting  the  initial  condition  x°  be  a  random  vector  and  then 
taking  the  expectation  of  the  relative  contribution  of  a  mode  to  a  state.  The  initial  condition 
is  assumed  to  be  a  random  vector  satisfying  the  following  assumption  of  independence. 

Assumption  1.  The  components  x®,  j  =  1  of  the  initial  condition  vector  x°  are 

independent  random  variables  with  probability  density  function  ■ 

3 

The  probabilistic  notion  of  participation  factors  for  participation  of  modes  in  states  fol¬ 
lows.  Denote  by  E  the  expectation  operator. 

Definition  3.  Suppose  Assumption  1  holds.  Define  p^,  the  participation  at  time  t  =  0  of 
the  mode  A,;  in  state  X/.,  as  the  expectation 


Pki  ■■=  E  { 


( llx0)r\ 


x 


o 

k 


} 


(15) 
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whenever  this  expectation  exists. 


In  applications  of  this  definition  to  specific  problems,  it  is  useful  to  rewrite  the  defining 
formula  (15)  as  follows: 


Pki 


E{ 


( llx°)rlk . 


Xi. 


ME 

3= 1 


n  (li~.0\„i 
\LjXj)'k 


X 


.0 


E{ 


T0  +  M  E 

3  =  1 


«+  E  biM-M 

i  =  i 

j  ^  k 


(16) 


To  verify  that  the  quantities  p*,;  of  Definition  3  do  not  depend  on  the  units  of  the  state 
variables  Xj,  j  =  1, ...  ,n,  invoke  Assumption  1  in  Eq.  (15)  and  expand  the  dot  product  in 
the  numerator,  to  get 


Pki 


^  1  ...o  I 


xi 


'  —oo  J  —  oo 


(liX3 


0  ^  nr>i 

•  u 


3  = 1 


XI 


)dx°1---dx°n_1dx°n.  (17) 


As  in  Section  2  (comments  following  Remark  1),  consider  a  change  of  variables  Xj  =  ocjXj, 
j  =  1, ...  ,n.  In  the  new  coordinates,  the  right  and  left  eigenvectors  corresponding  to 
eigenvalue  A*  have  components  r*  =  a3rl-  and  =  ay-1/*.  Also,  clearly,  the  initial  condition 
components  become  Xj  =  ajXp  j  =  1, . . .  ,n.  The  density  functions  for  the  initial  condition 
components  are  given  in  the  new  coordinates  by 

=  —  /*?(*?)•  (18) 

j  J  a.j  3 

The  differentials  scale  as 


daij  =  ajdxp  j  =  1, . . . ,  n. 


(19) 
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Denote  by  pki  the  participation  factors  computed  in  the  new  coordinates  using  Eq.  (17). 
Using  the  facts  noted  in  the  preceding  paragraph,  it  follows  immediately  that  pki  =  Pki ,  i.e., 
the  participation  factors  defined  in  Definition  2  are  indeed  independent  of  the  choice  of  units 
for  the  system  state  variables. 

The  following  well-known  fact  from  probability  theory  will  be  useful  in  the  examples 
below.  It  will  also  be  employed  in  Section  4  for  the  study  of  participation  of  states  in  modes. 

Lemma  1.  (See,  e.g.,  [2,  pp.  163-164].)  Let  X  and  Y  be  random  variables  and  let  g(X,Y ) 
be  a  function  of  X  and  Y.  Then 


E  {g(X,Y)}  =  EY{Ex{g(X,Y)\Y}  },  (20) 

where  the  notation  Ex  and  EY  is  used  to  emphasize  that  the  inner  expectation  is  conditioned 
on  Y  and  taken  with  respect  to  X,  and  the  outer  expectation  is  unconditional  and  taken 
with  respect  to  Y.  ■ 

For  the  purposes  of  this  paper,  the  following  observation,  based  on  Lemma  1,  is  particu¬ 
larly  valuable. 

Remark  3.  If  X  and  Y  are  independent  random  variables  and  the  probability  density  of  at 
least  one  of  X  or  Y  is  symmetric  with  respect  to  the  origin,  then: 


E  {XY}  =  0,  and 

(21) 

Eiy}  =  o- 

(22) 

Two  examples  of  the  use  of  Definition  3  are  given  next. 

Example  1.  Suppose  that  the  marginal  densities  fx°(xJ)  are  symmetric  with  respect  to 
x®  =  0,  i.e.,  that  they  are  even  functions  of  x°,  for  j  =  1, . . . ,  n.  With  this  assumption,  Eq. 
(16)  gives 


3  ±  k 


Jlrl 
Lk'  k 


(23) 
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by  Remark  3,  Eq.  (22).  Thus,  in  this  setting,  the  participation  factors  (15)  of  Definition  3 
reduce  to  those  of  the  original  definition  (3)  of  [9], [5], [10].  It  is  also  interesting  to  note  that 
the  calculations  above  are  still  valid  for  a  particular  state  Xk  for  which  the  corresponding 
initial  value  xk  is  known  to  be  distributed  symmetrically,  even  if  the  other  initial  conditions 
are  not  distributed  symmetrically.  ■ 

Example  2.  Suppose  that  the  state  variables  are  restricted  to  be  nonnegative,  and,  more 
specifically,  that  the  density  functions  fx0(x°j )  are  Rayleigh: 


< 


X0)2 

.7  ' 


if  Xj  >  0, 


0  if  Xj  <  0, 


(24) 


for  j  =  1 , ...  ,n.  Here  the  bj  are  positive  parameters  of  the  individual  Rayleigh  densities. 

The  mean  of  X®  is  E  { X j  —  1, . . . ,  n.  The  participation  factors  pki  of  Definition 
3  are  evaluated  using  (16)  as  follows  (Lemma  1  is  not  used  in  the  following  calculation, 
though  it  is  easy  to  check  that  it  would  give  the  same  result): 


'V* 
j'  k 


pH  =  iyk+  53  i: r 
j  =  i 
3  ±  k 


r°°  x°  x\x%  •  •  •  _(Al 

J  1  *  n  ^  9h-.  9 


'o  Jo  Jo  xl  6i&2  *  •  •  bn 


dx°x  ■  ■  ■  dx()n_1dx 


=  Ari+  E  lH 

3  =  1 

3  ±  k 


(‘oo  roc  • 


0\2  (*°)2  (*°)2 


*0  Jo  bjbk 


■e  i 


2b;  2  bt, 


k  dXjdxk 


Kr 


(xuy 
■j’ 

„0\2„  2b; 


n  poo  p  26,.  poo  (tv)2P 

4 r’  +  V  fir i  f  - dx«dx« 


°j'k 


3  =  1 
j  k 


=  44  +  E 

3  =  1 
j  ±  k 


(*l)2 

n  roo  p  2bfr 

Vi  [  — - E  {Xf}  dx°k 


j  k 


' 0 
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1  nn.1 

k 


44  +  E 

.7  =  1 

.7  ^  k 


I  nb. 


)dx°k 


7T  bi 


=  44  +  E  2 v/ 4 


.7  =  1 

.7  ^  k 


j  k ' 


(25) 


Thus,  the  participation  factors  obtained  for  this  example  differ  from  those  of  the  original  def¬ 
inition  of  [9], [5], [10].  This  illustrates  the  importance  of  initial  condition  uncertainty  models 
in  the  study  of  mode-state  interaction.  ■ 


4  Participation  Factors:  States  in  Modes 


In  the  foregoing  section,  the  emphasis  was  on  participation  of  modes  in  states.  In  the  present 
section,  attention  is  focused  on  assessing  the  other  aspect  of  mode-state  interaction,  namely 
participation  of  states  in  modes.  The  development  in  this  section  continues  to  address  the 
effect  of  initial  condition  uncertainty  on  mode-state  interaction.  The  uncertainty  is  modeled 
probabilistically  rather  than  set-theoretically,  since  the  authors  have  found  the  probabilistic 
approach  to  lead  to  simpler  analysis. 

Consider  again  Eq.  (1),  repeated  here  for  convenience: 


x(t )  =  Ax(t). 

Denote  by  V  the  matrix  of  right  eigenvectors  of  A: 

V  =  [r1  r2  ■  ■  ■  rn]. 


(26) 


(27) 


Recall  that  the  corresponding  left  eigenvectors  are  row  vectors  denoted  as  Z1,  l2, . . . ,  ln.  From 
the  normalization  (2)  and  the  orthogonality  of  left  and  right  eigenvectors  corresponding  to 
distinct  eigenvalues,  it  follows  that 


I/"1 


\lnJ 


(28) 
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Perform  the  change  of  variables 


z  :=  V  1x. 


(29) 


Then  z  follows  the  dynamics 


z(t)  =  V  1AVz(t) 
=  A  z, 


(30) 


where  A  :=  diag(Ai,  A2, . . . ,  An).  This  implies  that  the  new  states  zi:  i  =  1  evolve 

according  to 


Zi(t) 


=  eXitz° 


(31) 


where  z°  =  (z°, . . .  ,z°)T  :=  z(0).  It  is  clear  that  Zi(t)  represents  the  evolution  of  the  i-th 
mode.  To  dehne  the  participation  of  the  original  states  x k,  k  =  1, ...  ,n  in  zi:  z °  is  written 
in  terms  of  x°  and  then  an  appropriate  expectation  is  taken.  This  is  done  next. 

From  Eqs.  (29)  and  (28),  Zi(t)  is  given  by 

Zi(t)  =  eXitFx° 

n 

=  eA,‘E(tA)'  <32> 

3= 1 


This  equation,  which  shows  the  contribution  of  each  component  of  the  initial  state  j  = 
1  ,...,n,  to  the  i-th  mode,  motivates  the  following  definition  for  the  participation  factor 
governing  participation  of  states  in  modes. 

Definition  4.  Suppose  Assumption  1  holds.  Dehne  pki,  the  participation  at  time  t  =  0  of 
the  state  Xk  in  the  mode  A*,  as  the  expectation 

pH  :=  E{1^}  (33) 

Zi 

whenever  this  expectation  exists.  ■ 

Eq.  (33)  can  be  rewritten  in  a  form  that  can  be  directly  applied  depending  on  the 
probability  distribution  of  the  initial  condition,  much  in  the  same  way  as  was  done  for  Eq. 
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(15)  of  Definition  3.  To  wit,  rewrite  x k  in  terms  of  the  components  of  z°.  From  Eq.  (29),  it 
follows  that 


x°  =  Vz°. 


(34) 


Thus, 


x 


o 

k 


Erk?- 


3=1 


(35) 


Using  Eq.  (35)  in  Definition  4  gives  the  following  result,  analogous  to  Eq.  (16)  for  the  case 
of  participation  of  modes  in  states: 


Phi 


li  ™0 

rp  rLk±k  \ 

E  i  „0  J 


li  SZn  r3 
T  Z^j=i '  k.j 


=  E{ 


=  E{ 


Lk'  kzi 


}  + 


E  KeU} 

3  =  1 
3  ±  i 


1  (.  r[.  + 


.z? 


3  =  1 

j  #  « 


(36) 


The  general  result  (36)  along  with  Lemma  1  imply  that  Definition  4  reduces  to  the 
original  definition  of  participation  factors  (3)  under  Assumption  1  and  under  the  assumption 
of  symmetric  probability  densities  for  the  initial  condition  components  x°k.  This  is  illustrated 
in  the  following  example,  which  is  analogous  to  Example  1.  If  the  densities  for  the  initial 
condition  components  xk  are  not  symmetric,  then  the  corresponding  densities  for  the  z(- 
would  need  to  be  calculated  and  used  in  the  formula  (36). 

Example  3.  Suppose  that  the  marginal  densities  fxr>  (E-)  are  symmetric  with  respect  to 

3  J 

Xj  =  0,  i.e.,  that  they  are  even  functions  of  x°,  for  j  =  1 ,n.  Since  z°  is  given  by  a 
linear  transformation  of  x°,  E  {x0}  =  0  implies  that  E  {z(y}  =  0  also.  Therefore,  under  the 
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symmetry  assumption,  Eq.  (36)  along  with  Remark  3,  Eq.  (22)  give 


Pki 


lkrk  + 


n 


3  =  1 
3  i 


Ae 


nri 
"k  k- 


(37) 


5  Conclusions 

A  new  approach  to  the  dehnition  of  participation  factors  has  been  given.  The  approach 
does  not  assume  any  particular  choice  for  the  system  initial  condition.  Rather,  the  initial 
condition  is  taken  to  be  an  uncertain  quantity  and  the  participation  factor  is  defined  as  an 
average  contribution  of  a  mode  to  a  state,  or  of  a  state  to  a  mode.  The  averaging  is  taken  over 
the  set  of  possible  initial  conditions,  in  either  a  set-theoretic  or  a  probabilistic  framework. 
The  results  of  this  paper  were  shown  to  generalize  the  original  dehnition  of  participation 
factors  as  given  in  references  [9], [5], [10].  This  work  can  lead  to  further  useful  generalizations 
of  the  participation  factors  concept  in  various  directions.  For  instance,  generalizations  to 
other  system  types  can  be  pursued.  Also,  definitions  along  the  same  lines  could  facilitate 
analytical  treatment  of  the  relation  between  control  inputs  and  states/mode  interaction, 
which  in  turn  would  have  implications  for  actuator  and  sensor  placement.  The  authors  are 
pursuing  such  generalizations. 
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